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Abstract 

We present the results of the next-to-next-to-leading order QCD analysis of the recently 
revised experimental data of the CCFR collaboration for the xF^ structure function using 
the Jacobi polynomial expansion method. The effects of the higher twist contributions 
are included into the fits following the infrared renormalon motivated model. The special 
attention is paid to the checks of the predictive abilities of the infrared renormalon model 
and to the independent extractions of the x-dependent shape of the twist-4 contributions 
to the xF^ structure function in the process of the leading order, next-to-leading order 
and next-to-next-to-leading order fits of the revised CCFR data. We stress that at the 
next-to-next-to-leading order the results for as{Mz) turn out to be almost nonsensitive to 
the higher-twist terms. We obtain the following result {Mz) = 0.117±0.002(stat)± 

0.005(s?/st) ± 0.003(t/ieory). The comparison of the outcomes of our next-to-leading order 
and next-to-next-to-leading order analysis indicate that the theoretical QCD uncertainties 
were underestimated in the process of the next-to-leading order determination of as{Mz), 
made recently by the CCFR collaboration itself. 

Presented in part at the QCD session of Recontres des Moriond-97 (Mareh, 1997). 


1. The detailed QCD studies of the behavior of the structure functions of the non-polarized deep- 
inelastic scattering (DIS) are still very important both from the phenomenological and theoretical 
points of view (see e.g. the reviews ii). The revision of the experimental data of the CCFR 
collaboration, previously presented in Ref.Q put on the agenda the reconsideration of the results of 
the fits of the old CCFR data, performed at the next-to-next-to-leading order (NNLO) of perturbative 
QCD both without higher twist [|| and with higher twist contributions Q using the Jacobi polynomial 
expansion method. The results of application of this method for the extraction of the x-shape of the 
twist-4 contributions from the NNLO fits of the combined xFs data (besides the ones analyzed in 
Ref.Q) can be found in Ref.|^. The Jacobi polynomial approach was developed in Refs.[P|-[p^ 
and successfully used in the process of the next-to-leading order (NLO) hts |ll|- [ll4|| and NNLO hts 
[|I5| , ^ of different experimental data. 

R should be stressed that at the NLO level the results obtained turned out to be in good agreement 
with the outcomes of applications of the traditional variant of the DGLAP method 


solution of the integro-differential equation (compare the results of Refs.[]R 
neutrino-nucleon DIS and of Refs. 


and Ref. 


based on the 
in the case of 


12] with the ones of Ref. [jlTl ] in the case of the analysis of the data 
for the F2 structure function, obtained by the BCDMS and SLAC collaborations). However, the Jacobi 
polynomial expansion method is also giving the possibility to include into the QCD fits the existing 
information about the NNLO approximations of the coefficient functions and anomalous dimensions 
of the finite number of Mellin moments of the DIS structure functions, calculated analytically in 
Refs.|R]- 1^] and Ref.j^ respectively (for the detailed methodological description of the application 
of the Jacobi polynomial expansion machinery at the NNLO see Ref.[^). Therefore, in view of the 
current lack of information about the NNLO corrections to the Alarelli-Parisi kernel, the Jacobi 
polynomial expansion method turns out to be the unique method of including selfconsistently the 
available NNLO perturbative QCD effects into the fits of the DIS structure functions data. 

In the present work we demonstrate once more the compatibility of the DGLAP method and the 


( 4 ^ 

MS 


Jacobi polynomial expansion approach at the NLO level and extract the values of the parameter A 
from the revised CCFR data of Ref.Q (hereafter called as the CCFR’97 data) with taking into account 
target mass corrections and twist-4 effects, which were introduced in our analysis using the concept of 
the infrared renormalon (IRR) model of Ref. 122 j. At the next stage we include into our analysis the 
NNLO perturbative QCD corrections. We extract the x-shape of the twist-4 contributions to the xF^ 
structure function at the leading order (LO), NLO and NNLO and compare it with the IRR model 
predictions of Ref.jQ. We also get the values of as{Mz) both at the NLO and NNLO and compare 
them with the ones obtained by the CCFR collaboration in Ref.[^. 

2. The application of the Jacobi polynomial expansion machinery at the NLO and NNLO was 
described in detail in Refs.[^-[]R| and Refs.fT^, correspondingly. Nevertheless we will introduce 
again some basic definitions. The NNLO QCD Q^-evolution of the moments of the xF^ structure 
function = Jq x"' ^F^{x,Q'^)dx (where n = 2,3,4,...) is determined by the solution of the 

renormalization group equation as 
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Hn{Ql,Q^) ,n>2 


( 1 ) 


where 

[1 +p{n)As{Ql) + q{n){As{Ql)y] cJ5(yl,(Qg)) 


( 2 ) 


and C^]i!^l{As{Q‘^)) = I + C^^'>{n)As{Q'^) + C^‘^\n)A‘^^{Q‘^), As{Q‘^) = as{Q‘^)/{^7r), p{n), q{n) are 
related to the coefficients of the anomalous dimension functions 'j^g{As) and the QCD /3-function by 
the equations defined in Ref.[^ and the numerical expressions of the coefficients In's xFsi''^)^ 7.^5(^)’ 

(^) 5 (^) for / = 4 are given in our previous paper of Ref. j^j . 

Having at hand the QCD expressions for the Mellin moments M„(Q^) we can reconstruct the 













structure function xF^{x,Q‘^) using the Jacobi polynomial expansion method: 




^max ^ ^ 

=x"(l-x)^ £ 0"’^(x) ^ c5’"^(a,/3)M,+2 .xF3 (Q^) , 

n=0 


j=0 


( 3 ) 


where 0^’^ are the Jacobi polynomials and a, (3 are their parameters, fixed by the condition of the 
requirement of the minimization of the error of the reconstruction of the structure functions (see 
Ref.|W] for details) and Mn{Q‘^) is the QCD expression for the Mellin moments, corrected by the target 
mass contributions and twist-4 effects. Throughout this paper we will use the following expression 


MniQ'^) = + 


n{n + 1 ) Ml 


u -|- 2 Q 


+ c{n)Mi:^{Q^y-^^ + O(^) 
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where ^42 is the free parameter, M^ud is the mass of the nucleon, C{n) = —n — 4 -|- 4/n -|- 2/(n -|- 1) -|- 
4/(n + 2) + 4S'i(n) {Si{n) = Z)j=i 1/j) is determined by the IRR model calculations of Ref.@ (for 
the related studies of the IRR model predictions for the non-singlet structure functions of DIS see 
Ref. |]^ ). It is worth to mention here that the IRR model calculations are known to be complimentary 
to the ones, based on the dispersion approach for the definition of the running coupling, developed 
e.g. in the papers of Ref. [^1 . Note also, that the effects of the 0(l/(5^)-contributions in Eq.(4) come 
from two sources-the target mass effects (see Ref.^]) and the twist -6 contributions, which can be 
extracted from the IRR model calculations of Ref. |^] . However, the concrete fits of the experimental 
data demonstrate that the 0 (l/( 5 ^)-effects of the target mass corrections are negligibly small |5|. 
The similar feature was revealed while taking into account twist -6 contributions in the process of the 
fits of the CCFR’97 data (see the second citation in Ref.[^). Therefore, in our concrete studies we 
will neglect the 0(l/(5^)-contributions as the whole. We will also neglect the problem of defining 
selfconsistently the perturbative contribution to the coefficient function of the twist-4 correction in 
Eq.(4), which generally speaking differs from {Cp). 

One more theoretical input of our analysis is the factorization-scale dependent factor M,^‘^(Qg) in 
Eq.(l). Throughout this work we will use the MS'-factorization prescription, fixing [Ql) as 

dxx- 2 H(g 2 )/(Qg)(i _ ^)c(Qg)(i + (5) 

Jo 


where Qq is the relevant factorization scale. Note also that we will take in Eq.(3) Nmax = 10 for 
performing LO and NLO fits, while our NNLO fits will be done in the case of N^ax = 6 . At the NNLO 
we will use the calculated in Ref.|^] NNLO contributions 7 ^g( 2 ), 77 ^ 5 ( 4 ), 7 ^^( 6 ), 77 ^ 5 ( 8 ), 7 ^^( 10 ) 
to the anomalous dimension functions and the estimated by the procedure of smooth interpolation 
results for 7 ^ 5 ( 8 ), 7 ^ 5 ( 6 ), 7 ^ 5 ( 7 )) 7_^s(9) P- E is worth to mention here that the similar procedure 

was previously applied in the process of the NNLO non-singlet fits of the BCDMS data in Ref. [IT^ . 

f 2') 

However, we will avoid the application of the extrapolation of the expression for 7 )^ 5 (n,) at n > 10. 

To extract the x-shape of the twist-4 contributions we will simulate the xF^ structure function by 
the following expression 

h{x) 


xF^x, Q'y = xFy^{x, Q'y + 




( 6 ) 


where the Q^-dependence of the leading twist term xF^^(x, Q^) in Eq.( 6 ) is determined by perturba¬ 
tive QCD (including the proper treatment of the target mass corrections) and the problem of fixing the 
g^-dependence of the coefficient function of the twist-4 contribution is not considered. The constants 
h{xi) (one per each x-bin) are parametrizing the x-dependence of the twist-4 contribution, which in 
accordance with the theoretical predictions of Refs.[p5|, in the limit x —>■ 1 can be modeled as 
h{x) —> (1 — x)‘^“^, where the parameter c is introduced in Eq.(5) and like the parameters h{xi), A, b, 7 
and the QCD scale parameter A^^ can be determined by fitting the concrete sets of the experimental 
data for the xF^ structure function. Note, that the parameter A^^ is related to the expression for 
the QCD coupling constant Ug via the standard expansion in inverse powers of ln{Q’^/. 







2. At the first stage we performed the fits of the CCFR’97 data of Ref.[^ for the cuts of these data 
at a: < 0.7 using the value of the factorization point Qq = 5 GeV^ with the help of two independently 
written computer programs, which realize the Jacobi polynomial expansion method. The first Program 
has been created as the result of the works of Ref, [p^] , while the second Program previously found its 
more distinguished application in the process of the fits of Ref, [p!^] . The Jacobi polynomial parameters 
were fixed as a = 0.7 and /? = 3 like in our previous paper of Ref.Q. Both Programs are giving the 
identical results, which are presented in Table 1. All fits were done in the LO, NLO and NNLO 
level both without twist-4 contributions and with twist-4 corrections, taken into account through the 
expression of Eq.(4). 

In Table 2 the NLO results, obtained with taking into account the IRR-model estimates of the 
twist-4 contributions, are compared to the outcomes of the independent fits of the CCFR collaboration, 
which are presented in the second reference from Ref. Q. Let us remind that these results were got 
using the DGLAP-equation. One can see that within the statistical uncertainties our NLO results 
agree with the results obtained by the CCFR collaboration with the help of another method and 
another computer program. 

Looking carefully on Table 1 we arrive to the following conclusions; 

• The effects of the NNLO perturbative QCD contributions turn out to be very important in the 
analysis of the CCFR’97 data for the xFa structure function. Indeed, for the different Q^-cuts 
they are diminishing the values of the QCD scale parameter A^^ by over 50 — 130 MeV (!) 
provided the twist-4 corrections are taken into account through the IRR-model of Ref. ; 

• For the low Q^-cuts > \ GeV^ and > 5 GeV^ we observe the existence of the following 
hierarchy {x^)nnlo < (x^)iVLO < (x^)lO) which demonstrate the importance of taking into 
account of the effects of perturbative QCD corrections in the process of the fits of the concrete 
experimental data for the structure functions of DIS; 

• When the effects of higher-twist corrections are included, the NNLO results for A^^ are less 
sensitive to the values of Q^-cuts of the CCFR’97 experimental data than the similar NLO 
results; 

• The inclusion of the twist-4 contributions into the fits of the CCFR’97 data results in the 
detectable increase of the central values of A.^^ at the NLO for almost all Q^-cuts besides the 
cut > I GeV^ (to our point of view, the analysis of the data in the latest region can be 
more sensitive to the change of the value of the factorization scale Qq = 5 GeV^ say to the 
value Qq = 1 GeV^, which is more appropriate for fitting the data in the region of momentum 
transfered > 1 GeF^); 

• At the NLO level the values of the fitted parameter responsible for the inclusion of the 
twist-4 corrections, do not contradict the value A '2 ~ —0.2 GeV^, which was chosen in the work 
of Ref.p2| for the parametrization of the twist-4 corrections within the IRR model approach; 

• We found that at the NNLO level the IRR-motivated expression of the twist-4 contributions 

influence the outcomes of the fits less vividly than at the NLO level. This feature is related 
to the fact, that the central values of the parameter which result from the NNLO fits 

turn out to be rather small. Moreover, we even are unable to fix the sign of at the 

NNLO level, since the statistical error-bars of the value of this parameter are large (!!). This 
observed property might reveal the feature of the partial cancellation of the IRR contribution to 
the twist-4 correction of the DIS structure function with the possible part of the 0(l/(5^)-term, 
which can arise from the summation of the ultraviolet renormalon effects, discussed recently 
e.g. in Ref. p^ and partly taken into account in our analysis through the certain sets of the 
perturbative NNLO corrections. 


However, one should keep in mind the results of the recent quantum mechanical consideration of 
Ref.|^], which are demonstrating that the renormalons might give a noncomplete information on the 



structure of nonperturbative effects. In its turn this can indicate that the IRR higher-twist contribution 
to Eq.(4) might have non-controllable theoretical uncertainties. These interesting conclusions make 
very important the direct extraction of the information about the higher-twist contributions from the 
analysis of the CCFR’97 experimental data following the lines of the work of Ref.|p, devoted to the 
similar analysis of the old CCFR data. 

3. In this section we follow the ideas of Ref. 117, u and undertake the simultaneous extraction 
of the parameters A, b, c, 7 , and the x-shape of the twist-4 corrections from the CCFR’97 xF^ 
experimental data using the model of Fq.( 6 ). In accordance with the x-bin structure of the CCFR’97 
data we consider 16 bins Xj, presented in Table 3. We have taken the number of active flavours 
/ = 4 in all region of momentum transfered 5 GeV^ F < 199.5 GeV^. The results of the LO, 
NLO and NNLO fits, performed in the kinematical conditions Qq = 5 GeV^ and > 5 GeV^, are 
presented in Tables 3,4 and Figs.1-3. At Fig.2 we also plot the prediction for h(x), obtained within 
the framework of the IRR model approach of Ref.[22|. Note, that in the work of Ref. p^ the concrete 
predictions were obtained for the twist-4 contribution, written down in the slightly different form: 
xF^'^{x,Q‘^) = xF^"^D{ x)/Q‘^. To transform this prediction to our normalization of Fq.( 6 ) we used 
the program to compute the x-shape of this power term, written down by the authors of Ref.l^l^. 

Let us now describe the main conclusions, which follow from the results of Tables 3,4. 


• The statistical uncertainties of the values of A^^ are increasing drastically. This property 
results from increasing of the number of the fitted parameters and reflects the correlation of the 
uncertainties of the values of the parameters A, b, c, 7 and A^^ with the ones, which appear in 
view of the lack of the precise theoretical information about the twist-4 contributions; 


• The non-normalized to the number of degrees of freedoms value of the y^-parameter drastically 
decreases and turn out to be rather stable to the inclusion of the NLO and NNLO corrections 
into the fits of the CCFR’97 experimental data. This property is reflecting the positive feature 
of the incorporation of the additional 16 parameters h{xi) into the analysis of the CCFR’97 
experimental data; 


• We confirm the conclusion of Sec.2 about the importance of taking into account NNLO contribu¬ 
tions to Fq.(2) in the concrete fits of the experimental data. Indeed, the NNLO QCD corrections 
are decreasing both the central value of the parameter A^^ and the corresponding statistical 
uncertainties. Moreover, making additional fits, we found that the strong decrease of the value 
of the parameter A^^ at the NNLO is resulting from the incorporation of the NNLO corrections 

into Fq.(2), but not into the relation between q;s(Q^) and A.^^; 

IVl o 

• We have checked that the results of the LO and NLO fits are almost non-sensitive to the changes 
of Nmax = 10 in Fq.(3) to N^ax = 6 , considered in the process of the NNLO fits. This feature is 
rather welcomed from the point of view of the confirmation of the reliability of the application 
of the Jacobi polynomial expansion machinery; 


At the NLO and in the region of the intermediate values of x the form of the twist-4 correction 
h{x) turns out to be in satisfactory agreement with the x-behavior of the IRR-model prediction 
of Ref.|^]. This property is confirming the conclusion of Ref.pl 


that in the first approximation 
the IRR model does in fact good job as far as x-dependence of the higher-twist corrections is 
concerned; 


However, we also confirm the conclusion of Refs. 128, 291 that the IRR model might give the non- 
complete estimate of the mass scale of the higher-twist contributions. Indeed, we observe that 
the inclusion of the NNLO corrections into the fits of the experimental data is resulting in the less 
vivid x-dependence of the twist-4 correction h{x) (!), which model the effects of the truncation of 
the perturbative QCD expressions for the related moments and of the reconstructed 

structure function xF^^{x,Q‘^) at the definite order of perturbation theory. Moreover, looking 


'^We are grateful to B.R. Webber for providing us this program 











carefully on Table 3 we observe that in almost all Xj-bins the following tendency < 

{xi)\NLO is revealed, where h'^‘"{xi) denotes the central values of the corresponding twist-4 
parameters. These observations are one of the main outcomes of our analysis, presented in this 
Section. Following the hypothesis, made in the item 6 of Sec.2, we are relating the explanation 
of this feature to the possible NNLO manifestation of the summed up ultraviolet renormalon 
contributions, discussed in detail in Ref.|^]; 


• In the region of the intermediate values of x the results of our extarctions of h{xi) (see Figs.1-3) 
do not contradict to the predicted in Ref.[p^, behaviour |/i(x)| —> (1 — cY~^ (where we find 
3 < c < 4). Note, however, that the predictions of Ref.[ 
limit Xi 1, while our results are related to the kinematical region 0.0125 < Xi < 0.65, which 
in fact is lying quite far from this asymptotic regime. 


25, M were obtained in the asymptotic 


• The results of the fits demonstrate that the twist-4 contributions can have the sign-alternating 
character, discussed for the first time in the process of the fits of the DIS experimental data in 
Ref.[^], where another phenomenological approach was used for modeling higher-twist contri¬ 
butions. R should be stressed, that our picture of the x-behaviour of h{x) is more rich. Indeed, 
while the considerations of Ref, [p^] indicate the existence of only one zero in the values of h{x), 
we see that h{x) can have two zeros: the first one at small values of x and another one in the 
region of larger values of x (see Fig.1-Fig.3). 

4. In Ref.[^ the final result for the parameter extracted from the xF^ CCFR’97 data, 

was presented for the case of the fits, made within the kinematical conditions Qq = 5 GeV^ and 
> 5 GeV^. It was obtained with the help of the newly proposed non-conventional procedure of 
the treatment of the systematic uncertainties, called as the ’’global systematic ht”. It is based on the 
incorporation of the systematic variations of the structure functions as the part of the fit and on the 
redefinition of the x^-parameter. The application of this procedure results in the substantial reduction 
of the value of the newly defined x^-parameter and of the combined statistical and systematic error 
bars. However, in order to get the feeling, whether the theoretical error of the NLO value of was 
estimated in the process of the work of Ref.[^ correctly, we will compare the outcomes of our NLO 
and NNLO studies to the following CCFR’97 result, obtained in the second work from Ref.|^: 

NLO: A^ = 387 ±42{stat)±93{sy St) ± 17 {HT)±70{scale) MeV (7) 


where the statistical and systematic uncertainties were determined by the conventional way and the 
higher-twist error was estimated by varying the IRR model parameter A '2 in the following interval 
(0, —0.2 GeV"^) with the central value = A' 2 {DW )/2 = —0.1 GeV"^ ^ 

Notice, that in Eq.(7) the scale-scheme uncertainties, which in part reflect the sensitivity of the 
outcomes of the hts to the non-taken into account higher order QCD corrections, were ad hoc taken 
from the estimate Aas{Mz) = ±0.004, made in the process of the NLO studies of Ref.]^, aimed to 
the NLO analysis of the BCDMS F 2 structure function data of the fiN DIS. Evolving now the result of 
Eq.(7) through the threshold of the production of the 6-quark with the help of the derived in Ref.[31| 
relations with the choice rrn, ~ 4.6 GeV, we get the following value of as{Mz): 


NLO : as{Mz) = 0.122 ± 0m2{stat) ± 0m5{syst) ± O.OOl(HT) ± 0.004(sca/e) . (8) 


which corresponds to analysis of the CCFR’97 xF^ experimental data, performed in the second work 
of Ref.|^] with the separate determination of the statistical and systematic experimental uncertainties. 

Since we are able to incorporate in our fits the effects of the NNLO corrections and extract the 
value of A 2 - parameter (see Sec.2), we can get the feeling on the validity of the estimates of the 
theoretical uncertainties in Eq.(7) and Eq.(8). Indeed, for Qq = 5 GeV^ and > 5 GeV^ our NLO 

^The abbreviation A^DW) corresponds to the exact parametrization of the twist-4 contribution through the IRR- 
model of Ref. 







and NNLO results, which follow from the ones of Table 1, read 

NLO HT of Ref.[22] : = 371 ± 31{stat) ± 93{syst) MeV 

NNLO HT of Ref.r~ 


( 9 ) 


Ag_ = 293 ± 29{stat) ± 93{syst) MeV 


where the systematic uncertainties are taken from the original analysis of Ref.[^ (which we are unable 
neither to check nor to improve) and the statistical uncertainties in Eq.(9) are taking into account the 
deviation of the twist-4 parameter A '2 in Eq.(4) from its canonical value of Ref.[^|- 

One can see that the inclusion of the NNLO QCD corrections into the fits of the CCFR’97 ex¬ 
perimental data decreases the central value of the parameter by over 80 MeV. This amount is 
comparable with the estimate of the scale-scheme dependence uncertainty in Eq.(7). In its turn, this 
fact means that the theoretical uncertainties in the NLO results of Eq.(7) and Eq.(8) are underesti¬ 
mated. 

Indeed, it was already demonstrated in the second work from Ref.[^ that the outcomes of the 
NLO fits of the CCFR’97 experimental data for the xF^ structure function turn out to be sensitive to 
the upper value of the x-cut. The variation of this cut from x < 0.7 to x < 0.6 and to x < 0.5 (which 
is equivalent to neglecting the high Q^-data for the xEs structure function) is decreasing the central 
value of the CCFR’97 NLO result for A^^ by over 32 MeV and 40 MeV correspondingly^. Moreover, 
for the cut x < 0.7 there is still the place for other theoretical uncertainties (besides the controlled 
by us effects of the NNLO corrections), which appear due to the factorization scale dependence, the 
unknown N^LO contributions and the sensitivity to nuclear effects. Since the CCFR’97 data for the 
(xTa)^® structure function were not corrected by the latter ones, we made the fits, analogous to the 
ones described in Sec.3, but with modeling the ratio Rpe = by the ratio Rd = F^/F^, 

which was fixed by the approach, considered in Ref.|^]. These fits also result in the effect of the 
decreasing of the value of the parameter A^^, but by less considerable amount of over 15 MeV. It is 
highly desirable to check this conclusion by using some other approach for modeling the value of Rpe^ 
say the one, developed in Ref.^^ in the case of the analysis of the F 2 structure functions data. 

Evolving now the results of Eq.(9) through the thresholds of the production of the 6-quark and 

/ c'\ 

substituting the values of the parameter ALM into the inverse-log expression for as{Mz), we get 

IVl o 


NLO HT of Ref.[22] : 
NNLO HT of Ref. I 


as{Mz) = 0.121 ± 0.002(stat) ± 0.005(syst) ± 0.006(t/ieory) (10) 
as{Mz) = 0.117 ±0.002{stat)±0.005isySt) ±0.003{theory) . 

The theoretical error bars of the NLO and NNLO values of as{Mz) in Eq.(lO) are taking into account 
all discussed uncertainties plus the arbitrariness of application of the procedure of passing threshold 
of the production of 6-quark, which following the considerations of Ref.|M] we estimate as Aas{Mz) ± 
0 . 001 . 

The NNLO results from Eq.(lO) turns out to be in agreement with the following NNLO value of 
as{Mz). as{Mz) = 0.115 ± O.OOl(stat) ± 0.005(syst) ± 0.003{twist) ± 0.0005(sc/ieme), which was 
extracted in Ref. [^] from the comparison of the old CCFR data for the Gross-Llewellyn Smith sum 
rule at the reference point = 3 GeV^ |35| with the results of the perturbative NLO and NNLO 
calculations of Ref.|36|, supplemented with the information about the concrete values of the twist-4 
corrections, obtained in Ref.[^]. The optimistic scheme-dependence error in the NNLO result of 
Ref. |]^ is reflecting the uncertainty due to the unknown effects of the higher N^LO perturbative 
QCD corrections, estimated later on in the works of Ref. [08| . It should be stressed, however, that in 
view of the changes of the CCFR uN data the experimental analysis of Ref. [^] should be updated. 

Let us also mention, that the discussed in Sec.3 fits of the CCFR’97 data, which have the aim to 
extract the x-shape of the twist-4 contributions and to verify the predictions of the IRR-model, give 
as the outcomes the comparable values of A.^^ and thus as{Mz). Indeed, when transformed through 
the thresholds of the production of the 6-quark, the NLO and NNLO resnlts of Table 4 


NLO HT free : 
NNLO HT free : 


= 428 ± 158(exp) MeV 

IVl o 

A^^ = 264 ± 85(exp) MeR 

IVl o 


( 11 ) 


^Note, that for the cut x < 0.6 we confirmed this decreasing tendency. 








are equivalent to the following values of as{Mz)- 


NLO HT free: as{Mz) = 0.12A ±0.007(exp) ±0.0W{theory) (12) 

NNLO HT free : as{Mz) = 0.115 ± 0.006{exp) ± 0.003{theory) , 


where the experimental uncertainties are estimated by the statistical errors (which are increased due 
to the inclusion of the additional 16 twist-4 parameters h(xi) into the fits), the NLO theoretical errors 
are fixed by the difference between the NLO and NNLO results and the remaining NNLO theoretical 
ambiguities are determined by the uncertainties, taken by us into account in the process of estimates 
of the theoretical ambuguities of the NNLO value of as{Mz), presented in Eq.(lO). 

It should be stressed, that our results of Eqs.(10),(12) do not contradict to the NLO value as{Mz) = 
0.113 ± 0.003(ea;p) ± O.OOi^theory), obtained in the process of the NLO analysis of the BCDMS and 
SLAG data for the F 2 structure function of pN and eN DIS |17|. However, there are still the number of 
questions to be answered in future. The first, theoretical one, is related to the necessity of the inclusion 
into the analysis of the BCDMS and SLAG data the NNLO perturbative QCD corrections to the singlet 
part of the F 2 structure function (we consider the results of the NNLO non-singlet fits of Ref.Jl^ still 
non-complete, since in the process of these fits only the data for x > 0.35 was analyzed without 
taking into account target mass corrections). The second, pure experimental, question is related to 
the problem whether the scientific community will accept the proposed in Ref.|p new approach to 
treating experimental uncertainties. Note, that the accepted at present conventional way of their 
definitions is leading to the conclusion, that the combination of the traditionally defined statistical 
and systematic uncertainties of the GGFR’97 experimental data is not too small. Keeping this in 
mind, we agree with the conclusions of the analysis of the experimental data for the Bjorken polarized 
sum rule [39, 401 that in order to extract the more precise value of as{Mz) from the characteristics 
of DIS with taking into account both the higher-order perturbative QGD corrections and higher-twist 
effects it is necessary to have more precise experimental data. This conclusion is supporting the 
detailed investigation of the third problem, revealed in the process of the analysis of the GGFR’97 
experimental data of Ref.[^. It is related to the necessity of more detailed investigation of the possible 
discrepancies between the GGFR’97 F 2 uN data in the region of small x {x ^ 0.01 — 0.07) and the 
similar experimental data for the F 2 structure functions of the DIS of charged leptons on nucleons, 
measured by the BGDMS, NMG and SLAG collaborations. Among the possible ways of attacking 
this experimental problem are the possible more detailed studies of vAl DIS data of the IHEP-JINR 
neutrino detector (for the already available preliminary results see Ref.|41|) and the continuation of 
the careful measurements of the structure functions of vFe DIS at Fermilab Tevatron, which is the 
main goal of the NuTeV collaboration. 
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AgL (MeV) 

A 

b 

c 

7 

A'(HT) 

xVnep 

1 

LO 

277±31 

4.65±0.14 

0.69±0.01 

3.90±0.04 

1.82±0.13 


146.0/110 


367±43 

3.95±0.40 

0.63±0.03 

3.87±0.05 

2.48±0.47 

-0.25±0.04 

114.4/110 

NLO 

322±27 

3.67±0.72 

0.62±0.05 

3.76±0.10 

2.32±0.94 

- 

114.1/110 


322±23 

3.61±0.11 

o.eiio.oi 

3.76±0.04 

2.39±0.12 

-0.07±0.04 

111.4/110 

NNLO 

275±24 

4.01±0.34 

0.64±0.02 

3.61±0.06 

1.56±0.33 

- 

104.9/110 


276±25 

4.00±0.34 

0.64±0.02 

3.61±0.06 

1.56±0.32 

0.01±0.04 

104.9/110 

5 

LO 

266±37 

5.13±0.46 

0.72±0.03 

3.87±0.05 

1.42±0.33 


113.2/86 


436±56 

4.73±0.39 

0.68±0.03 

3.79±0.05 

1.82±0.31 

-0.33±0.06 

82.8^6 

NLO 

341±41 

4.05±0.38 

0.65±0.03 

3.71±0.06 

1.96±0.36 

- 

87.1/86 


371±31 

3.94±0.12 

0.64±0.01 

3.69±0.04 

2.05±0.12 

-0.12±0.05 

82.4/86 

NNLO 

293±29 

4.25±0.38 

0.66±0.03 

3.56±0.07 

1.33±0.33 

- 

78.4/86 


293±29 

4.25±0.37 

0.66±0.03 

3.56±0.07 

1.33±0.32 

-0.01±0.05 

78.4/86 

10 

LO 

287±37 

6.52±0.20 

0.79±0.01 

3.73±0.04 

o.eoio.io 


77.7/63 


531±74 

6.48±0.65 

0.79±0.04 

3.68±0.05 

0.94±0.26 

-0.52±0.11 

58.0/63 

NLO 

350±36 

5.29±0.21 

0.74±0.01 

3.64±0.06 

1.06±0.17 

- 

64.4/63 


439±55 

4.92±1.06 

0.73±0.07 

3.61±0.08 

1.42±0.75 

-0.24±0.10 

58.6/63 

NNLO 

308±34 

4.65±0.53 

0.70±0.04 

3.52±0.08 

1.06±0.38 

- 

58.7/63 


313±37 

4.64±0.52 

0.70±0.04 

3.52±0.08 

1.07±0.38 

-0.03±0.09 

58.6/63 

15 

LO 

319±41 

7.58±0.24 

0.84±0.01 

3.61±0.05 

0.16±0.09 


58.8/50 


531±57 

7.76±0.26 

0.86±0.02 

3.61±0.04 

0.45±0.09 

-0.57±0.13 

50.2/50 

NLO 

365±40 

6.28±0.21 

0.80±0.02 

3.60±0.04 

0.58±0.09 

- 

53.0/50 


441±44 

5.98±0.19 

0.80±0.02 

3.57±0.05 

0.87±0.10 

-0.25±0.13 

50.9/50 

NNLO 

314±37 

4.68±0.75 

0.70±0.05 

3.53±0.08 

1.08±0.51 

- 

50.9/50 


308±45 

4.69±0.72 

0.70±0.05 

3.53±0.08 

1.07±0.50 

0.03±0.14 

50.8/50 


Table 1. The results of the fits of the xF^ (i/Fe) revised data of the CCFR collaboration. The values 
of and A 2 {HT) are given in GeV^. The statistical errors are presented, {x^ is normalized to the 
number of experimental points (nep).) 


cut (GeV^) 

X cut 

Order 

(MeV) 

xVnep 

> 5 

X < 0.7 

NLO 

371 ± 30 

82.4/86 


CCFR result 


387 ± 42 

81.8/86 

g2 > 10 

X < 0.7 

NLO 

439 ± 55 

58.6/63 


CCFR result 


422 ± 47 

58.0/63 

g2 > 15 

X < 0.7 

NLO 

441 ± 44 

50.9/50 


CCFR result 


433 ± 55 

50.0/50 


Table 2. The comparison of the results of our NLO fits with the revised results provided by the 
CCFR collaboration in Ref. [3] 













LO 

NLO 

NNLO 

xVnep 

66 . 2/86 

65.6/86 

65.7/86 

A 

5.44 ± 1.74 

3.70 ± 1.56 

4.54 ± 0.88 

b 

0.74 ± 0.10 

0.66 ± 0.11 

0.69 ± 0.06 

c 

4.00 ± 0.18 

3.78 ± 0.21 

3.72 ± 0.19 

7 

1.72 ± 1.25 

2.86 ± 1.72 

1.43 ± 0.69 

Xi 

h{xi) \GeY^\ 

0.0125 

0.206 ± 0.321 

0.213 ± 0.332 

0.170 ± 0.302 

0.0175 

0.061 ± 0.268 

0.091 ± 0.289 

0.030 ± 0.241 

0.025 

0.146 ± 0.204 

0.220 ± 0.241 

0.136 ± 0.178 

0.035 

-0.021 ± 0.185 

0.114 ± 0.240 

0.008 ± 0.179 

0.050 

0.031 ± 0.142 

0.245 ± 0.230 

0.124 ± 0.164 

0.070 

-0.145 ± 0.127 

0.138 ± 0.233 

0.033 ± 0.165 

0.090 

-0.177 ± 0.125 

0.139 ± 0.225 

0.076 ± 0.166 

0.110 

-0.340 ± 0.126 

-0.015 ± 0.205 

-0.026 ± 0.166 

0.140 

-0.404 ± 0.114 

-0.092 ± 0.147 

-0.027 ± 0.141 

0.180 

-0.350 ± 0.164 

-0.077 ± 0.122 

0.054 ± 0.127 

0.225 

-0.554 ± 0.237 

-0.348 ± 0.170 

-0.167 ± 0.135 

0.275 

-0.563 ± 0.334 

-0.462 ± 0.272 

-0.196 ± 0.193 

0.350 

-0.314 ± 0.418 

-0.368 ± 0.371 

0.070 ± 0.200 

0.450 

-0.117 ± 0.415 

-0.266 ± 0.401 

0.183 ± 0.213 

0.550 

0.087 ± 0.333 

-0.109 ± 0.352 

0.097 ± 0.234 

0.650 

0.377 ± 0.215 

0.221 ± 0.244 

0.259 ± 0.188 


Table 3. The results of the LO, NLO and NNLO QCD fits of the CCFR’97 xF^ structure functions 
data for the values of twist-4 contributions h{x) and for the parameters A, b, c, 7 with the corresponding 
statistical errors. 



h{x) / 0 free 

h{x) ^ 0 of Ref. [22] 

h{x) = 0 


X^/nep 

{MeV) 

/nep 

A^ {MeV) 

^ Imp 


LO 

NLO 

NNLO 

66 . 2/86 

65.6/86 

65.7/86 

338 ± 169 

428 ± 158 

264 ± 85 

82.8/86 

82.4/86 

78.4/86 

436 ± 56 
371 ± 31 

293 ± 29 

113.3/86 

87.1/86 

78.4/86 

266 ± 37 

341 ± 41 

293 ± 29 


Table 4. The comparison of the results of the LO, NLO and NNLO fits, of the CCFR’97 xF^ 
structure functions data, performed with different assumptions of the values of the twist-4 parameter 
h{x). 










Figure captions. 

Fig-1 The results of the LO extraction of the x-shape of the twist-4 contribution h{x). 

Fig.2 The results of the NLO extraction of the x-shape of the twist-4 contribution h{x). For the 
comparison, the IRR-model prediction of Ref. [22], obtained using the NLO MRS parametrization, is 
also depicted. 

Fig.3 The results of the NNLO extraction of the x-shape of the twist-4 contribution h{x). 
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